www.gusucode.com > Weighted Differential Evolution Algorithm (WDE) > Weighted Differential Evolution Algorithm (WDE)/cc_wde/PointCloudProcessing/my_tin_simplifier.m

    function [out,fnc,v]=my_tin_simplifier(X,mydata)

n0=mydata.n;
xtest=mydata.xs;
ytest=mydata.ys;
ztest=mydata.zs;
fnctest=mydata.fnc;
N=size(X,1);
out=rand(N,1); % memory
for i=1:N
    
    
sol=X(i,:);
sol=reshape(sol,n0,2); % planimetric locations of tie points
zsol=fnctest(sol(:,1),sol(:,2));
v=[sol zsol];
fnc = scatteredInterpolant(v(:,1),v(:,2),v(:,3),'natural');
z0=fnc(xtest,ytest);
err=sum(abs(z0-ztest));

out(i)=err;


end